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Dynamics of coarsening of a statistically homogeneous 
fractal cluster, created by morphological instability of growth 
during an early stage of a first order phase transition, is in- 
vestigated theoretically. An exact mathematical setting of 
the problem, obeying a global conservation law, is presented. 
A statistical mean field theory is developed that accounts for 
shadowing and assumes that the fractal dimension of the clus- 
ter is invariant in time. The coarsening dynamics are found 
to be self-similar, and the corresponding dynamic scaling ex- 
ponents are calculated for any Euclidean dimension. 

PACS numbers: 64.60.-i, 6f .43.Hv, 64.60.Ak, 05.70.Fh 



A non-equilibrium dissipative system undergoes relax- 
ation to equilibrium after the driving agent is "switched 
off' or depleted. In complex systems the relaxation dy- 
namics can be quite complicated, and it is natural to 
seek for the dynamic scaling and universality. An in- 
structive, exactly solvable example of dynamic scaling in 
relaxation (coarsening) of rough surfaces with nonconser- 
vative dynamics is given by the deterministic (undriven) 
KPZ-equation |^,|^ . A much older example is decay of ho- 
mogeneous and isotropic hydrodynamic turbulence . 
Finally, there is an important class of relaxation prob- 
lems related to phase ordering dynamics, non-conserved 
and conserved, in the bulk [^[-Q and on the surface 

In the case of conserved dynamics, "switching off' 
of the driving agent occurs naturally as a result of a 
global conservation law. There is an important class of 
non-equilibrium systems that exhibit strong morpholog- 
ical instabilities and ramified growth at an early stage 
of the dynamics, show phase ordering at an intermedi- 
ate stage, and finally approach a simple thermodynamic 
equilibrium. Examples are provided by realistic diffusion- 
controlled growth systems, such as deposition of solute 
from a supersaturated solution and solidification from an 
overcooled liquid. The stage of morphological instability 
and its implications has been under extensive investiga- 
tion [p[-p^. On the contrary, the (usually much longer) 
phase-ordering stage has not been addressed. We claim 
in this Letter that this generic phase-ordering stage ex- 
hibits intrinsic self-similarity, calculate the scaling expo- 
nents for any Euclidean dimension d > 1 and make pre- 
dictions that can be tested in experiment or numerical 
simulations. 

We will concentrate, for concreteness, on the problem 
of growth of a single nucleus (we will call it cluster) from 
a supersaturated solution and start with a minimalistic 



mathematical setting that describes correctly the whole 
dynamics, from the stage of growth to the final equilib- 
rium. Let u(r, t) be the mass concentration of the so- 
lution normalized to the (constant) density of solute in 
the compact solid phase. The field u is governed by the 
diffusion equation 



du 

at 



(1) 



in a finite d-dimcnsional domain. We specify a no-flux 
boundary condition. 







(2) 



on the external boundary F, where index n stands for 
the normal component of a vector. Assuming that the 
moving interphase boundary 7 is in local thermodynamic 
equilibrium, we employ the Gibbs-Thomson relation 



Uq{1 + Xqk) 



(3) 



where uq is the (normalized) equilibrium concentration 
of the solution in the bulk, Xq is the capillary length and 
K is the local curvature for d — 2, or the mean curvature 
for d > 2. (For simplicity, we will limit ourselves to an 
isotropic surface tension.) Finally, mass conservation at 
the moving boundary yields the well-known relation for 
the normal speed: 



1 



u 
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It is easy to check that Eqs. (0)-(0) preserve the total 
mass of the solute. In the normalized form 



fir + I udr = const , 
'n 



(5) 



where flc is the volume (area) of the cluster, while 
denotes the region unoccupied by the cluster. This sim- 
ple conservation law describes the depletion of solution 
during the deposition process and imposes an important 
constraint on the dynamics. This constraint docs not ap- 
pear in the more traditional formulations of the diffusion- 
controlled growth problem ||9|-[r^, where an infinite sys- 
tem is studied, and the boundary condition correspond- 
ing to a constant (positive) flux or constant supersatura- 
tion at r ^ 00 is used instead of Eq. (||). Notice that, 
even in the limit of strong diffusion, it is the full diffusion 
equation (rather than its Laplace's equation limit) and 
no-flux condition on F that provide the conservation law. 
Also, the usually small term u in the denominator of Eq. 
(|) should be kept to get Eq. (|) right. 
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Despite a simple formulation, the moving boundary 
problem (0)-(|^) is formidable. In the remainder of this 
Letter we will present a mean field theory of coarsening 
that exploits a strong resemblance of this problem to the 
classical problem of Ostwald ripening [^-Q , accounts for 
shadowing and is based on the assumption that the frac- 
tal dimension of the coarsening cluster remains constant 
on an interval of scales shrinking with time [ p4|pq] . 

In a sufficiently large and "noisy" system P)-® with 
small initial supersaturation, a fundamental morpholog- 
ical instability ||l^ leads to ramification of the moving 
interface, so that the interface becomes fractal (DLA- 
like) . Because of mass conservation, the supersatura- 
tion decreases with time. Correspondingly, the surface- 
tension effects mitigate and finally switch off the mor- 
phological instability. The coarsening stage is character- 
ized by the (properly scaled down) supersaturation be- 
coming comparable to the mean curvature of the inter- 
face. The total cluster mass asymptotically approaches 
a constant value Mq , and the main effect at this stage is 
the diffusion-controlled mass transfer between different 
branches of the cluster that, in view of Eq. (||), pro- 
motes growth of larger branches (that is, those with a 
smaller curvature) at the expense of smaller ones (those 
with a larger curvature). This dynamics strongly resem- 
bles Ostwald ripening (OR) The classical problem 
of OR deals with statistics of many separate "drops" of 
the new phase, and the corresponding statistical mean- 
field theory assumes a negligibly small volume (or 
area) fraction of the drops, and completely ignores spa- 
tial correlations. It is very important that the dynamic 
scaling exponent predicted by the mean-field theory p^ ] 
agrees very well with experiments and simulations per- 
formed under conditions of a finite, or even large volume 
(area) fraction ^,||,D and significant spatial correlations 
(even when description in terms of individual drops be- 
comes irrelevant [^|). Inspired by this success, we propose 
a mean-field theory of coarsening for a statistically ho- 
mogeneous fractal cluster (that is a cluster whose fractal 
dimension is constant in space). 

Assuming for simplicity, that the cluster is (statisti- 
cally) azimuthally (for d — 2), or spherically (for d — 3) 
symmetric, we introduce a time-dependent distribution 
function f{R,r,t) for the number of branches with size 
R situated at the radius r from the cluster center. Look- 
ing for self-similarity, we assume a single scaling in R 
(which implies, in particular, that the typical branch 
size scales like its curvature radius) and a single scal- 
ing in r for the whole cluster. An additional vari- 
able is an effective (coarse-grained) supersaturation field 
Au{r,t) — u{r,t) — uq. The r-dependence in the fmic- 
tions / and Au accounts for shadowing and, as we will 
see later, for shrinkage of the coarsening cluster. This de- 
pendence has no analog in the "conventional" mean-field 
theory of OR @. 

Neglecting direct coalescence of branches, we write 



down a continuity equation 

f + |^(^«/) = 0, (6) 
while for Vr we adopt the well-known relation p7| , p^ 

following from the solution of an idealized problem of the 
Laplacian growth of an isolated spherical particle with 
the boundary condition (||) on its interface (all dimen- 
sional coefficients are scaled down) . Actually, the precise 
form of Eq. (|7|) is unessential for our purposes. What is 
really important (and used in the following) is that the 
characteristic normal speed of the boundary (that is, ex- 
pansion/contraction rate of an individual branch) scales 
like R~^, while the effective supersaturation scales like 
the characteristic mean curvature (see below) . 

Introduce the (time-dependent) mass content of the 
cluster within the radius r: 



M{r, t) 



dn'^ 



[r(| + 1)]^ 



nr poo 

■ / drr'^-^ / dRR'^ f{R,r,t). 
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(8) 



where r(a;) is the gamma-function. Mass conserva- 
tion of the coarsening cluster implies that M{rc,t) = 
Mq =const, where Tc is the (time-dependent) cluster ra- 
dius. 

We are looking for self-similarity: 

f{R, r, t) = t-^'^Rt-", rt'^) , Aw(r, t) = t-°'U{rt'^). (9) 

Correspondingly, the cluster radius will scale like rc(i) = 
riot~^, where 770 is a constant depending on d. Notice, 
that the anzats for / assumes self-similarity with respect 
to each of the two arguments, R and r. Upon substitu- 
tion, Eqs. (|) and (g) yield 



a = ly = 1/i 



(10) 



independently of the Euclidean dimension d and of (still 
undetermined) (3 and fj,. Employing the mass conserva- 
tion, we obtain a relation between f3 and /z: 



d+1 

dP + fi -— = . 



(11) 



Now we will use the assumption of invariance of the frac- 
tal dimension df [on an interval of scales [lmin(t)j Imaxit)] 
that is shrinking with time] . It is convenient to work with 
the mass density of the cluster: 



r(f + 1) 9M(r,i) 



dr 



r(i + i) 



(12) 
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where r] = rt^ . 

For a (statistically self-similar) "mass fractal" with the 
fractal dimension df one has M{r,t) — A{t)r'^f which 
implies a density 



p{r, t) ^ A{t) f 



df —d 



A{t) tP(d-df)^df-d ^ 



(13) 



[Since the coefficients of the scaling laws (and even the 
shape functions $ and U) will remain undetermined, we 
have omitted here and in the following all d- and df - 
dependent coefhcients.] A direct comparison of Eqs. (|l2) 
and ( p^ ) yields 



A(t) ^ and 

poo 

/ d^C'$(e,^)^'7'^"'- 
Jo 



(14) 



At scales smaller than Imin (t) the cluster is already com- 
pact. It is natural to identify Iminit) with a typical (for 
example, mean) value of the branch size R(t) that, ac- 
cording to Eq. (p^, scales with time like t^^^ There- 
fore, at radii r < R{t) the mass density must be constant. 
On the other hand, Eqs. ( p^ ) and ( p^ predict the fol- 
lowing behavior of the mass density at r = R{t): 



p{R,t) r^t'^~^'{Rt^ff 



^-{d-df)f3~tj.+ - 



Requiring this to be constant, we have 



{d-df)f3 + p 



df + 1 



(15) 



Combining Eqs. (^l]) and (|T^) we find 



^_ d-df _2d+l d 
3df ^ ^ 3 3df 



One can see that /3 > for any d, so the coarsening clus- 
ter is indeed shrinking with time . Also, /x > for any 
d. Having found the dynamic scaling exponents a, /?, fj, 
and V, we can make a number of predictions. A "global" 
prediction concerns scaling of the cluster perimeter (for 
d — 2) or interface area (for d — 3). These quantities de- 
crease like (independently of d and djl). Another 
"global" prediction concerns the already mentioned scal- 
ing of the cluster radius: Tc ~ t^id-df)/3df ^ p^j, example, 
for a DLA-like cluster {df ~ 1.71 for d ~ 2, and df ~ 2.5 
for d = 3), our predictions are Vc ^ ^-0.057 g^j^j ^-o.o67^ 
respectively. Returning to Eq. (|l^), we see that the 
mass density p{r,t) ~ t'^r'''f~'^, where a = {d — d/)/3. 
For d — 2 and 3 we get a ~ 0.10 and 0.17, respectively. 

Self-similar phase-ordering dynamics of a fractal clus- 
ter that we have considered crosses over to a compact re- 
laxation dynamics on a time scale for which Iminit*) ~ 
Imaxit*), that is R{t*) rc{t*). Finally, the cluster ac- 
quires spherical (circular) shape, corresponding to the 
only possible two-phase equilibrium in the system 12^ . 



Our mean-field theory is based on the assumption of a 
statistically homogeneous cluster. It is well known that 
the diffusion-controlled fractal growth crosses over to a 
compact growth when the growing cluster reaches the" 
diffusion length" Our estimates show that, for a suf- 
ficiently large system and small initial supersaturation, 
the depletion effects become important before the clus- 
ter radius has a chance to reach the diffusion length. As 
the result, the cluster stops growing and the coarsening 
effects that we have described starts dominating. 

A number of refinements of the theory are straightfor- 
ward (like an account for growth anisotropy by introduc- 
ing an additional, angle-dependent factor in the mean- 
field variables / and Aw). What is probably more impor- 
tant is that similar mean-field theories can be developed 
for other transport mechanisms (for example, interface- 
controlled coarsening, see, e.g. [^). The main difference 
will be in the value of exponent v, and scaling arguments 
yield v = 1/2 in this limit. In addition, we must put 
U{ri) — const, as the supersaturation field is uniform in 
this case . Once ly is found, the rest of exponents can 
be easily calculated. For the general v, one obtains 



^jd-df) 
df 



, , d-' 
v\2d+l~ — 

df 



Correspondingly, the interface area (perimeter in 2d) will 
decrease like t~'^ independently of d and df. 

Predictions of the mean-field theory can be checked 
in experiment and numerical simulations. Simulations 
are possible with different types of models. For the 
diffusion-controlled coarsening, one option is the sharp- 
interface model (|l|)-(3). Another can employ a conserved 
phase-field model As noise is not important at 

the coarsening stage, one can use the noiseless version 
of any of these two models and start directly from a 
fractal cluster [^. One more alternative would be to 
employ a modified DLA algorithm with the surface ten- 
sion psj] and total mass conservation properly taken into 
account. The simulations can be time-consuming, as the 
(pronounced) coarsening stage is necessarily much longer 
than the growth stage. 
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